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O ABSTRACT 

\^ The feedback between massive stars and the interstellar medium is one of the most 

^ important processes in the evolution of dwarf galaxies. This interaction results in nu- 

O merous neutral and ionised gas structures that have been found both in the disc and 

i in the halo of these galaxies. However, their origin and fate are still poorly understood. 

We here present new Hi and optical data of two Magellanic irregular dwarf galaxies 
^ in the Local Volume: IC 4662 and NGC 5408. The H I line data were obtained with the 

Australia Telescope Compact Array and are part of the "Local Volume H I Survey" . 
i_i They are complemented by optical images and spectroscopic data obtained with the 

ESO New Technology Telescope and the ESO 3.6m telescope. Our main aim is to 
study the kinematics of the neutral and ionised gas components in order to search 
for outflowing gas structures and to make predictions about their fate. Therefore, we 
perform a Gaussian decomposition of the H I and Ha line profiles. 
We find the H I gas envelopes of IC 4662 and NGC 5408 to extend well beyond the op- 
tical discs, with H I to optical diameter ratios above four. The optical disc is embedded 
into the central Hi maximum in both galaxies. However, higher resolution Hi maps 
show that the Hi intensity peaks are typically offset from the prominent Hll regions. 
While NGC 5408 shows a fairly regular H I velocity field, which allows us to derive a 
rotation curve, IC 4662 reveals a rather twisted H I velocity field, possibly caused by a 
recent merger event. We detect outflows with velocities between 20 and 60kms _1 in 
our Ha spectra of both galaxies, sometimes with H I counterparts of similar velocity. 
We suggest the existence of expanding superbubbles, especially in NGC 5408. This is 
also supported by the detection of FWHMs as high as 70kms _1 in Ha, which cannot 
be explained by thermal broadening alone. In case of NGC 5408, we compare our re- 
sults with the escape velocity of the galaxy, which shows that the measured expansion 
velocities are in all cases too low to allow the gas to escape from the gravitational 
potential of NGC 5408. This result is consistent with studies of other dwarf galaxies. 
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1 INTRODUCTION 

The interplay between massive stars and the interstellar 
medium (ISM) has a large effect on the formation and the 
evolution of galaxies, especially of dwarf galaxies. Photoioni- 
sation is the most likely interaction process. However, stellar 
winds and supernovae (SNe) explosions of massive stars also 
contribute significantly to the energy input into the ISM. 
The radiative and mechanical feedback heat the gas and 
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drive it outwards, sweeping up the ambient gas into a thin 
shell. A superbubble filled with hot gas evolves and begins 
to expand into the ISM. Due to Ray leigh- Taylor instabili- 
ties, the outer shell can rupture and the hot gas can vent 



out through so-called chimneys into the halo (|Norman & 
Ikeuchi||1989 ) 



Numerous ionised gas structures such as supergiant 
shells or filaments close to the galactic disc of dwarf galaxies, 
but also at kpc-distances from any place of current star for- 
mation, have been detected on deep Ha images (e.g., [Hunter 



et al.||1993| |Bomans et al |[l997l [Hunter fe Gallagher||1997 



Lopez-Sanchez &; Esteban 2008). Spectroscopic observations 
of the Ha line revealed that most of the ionised gas struc- 
tures expand from the disc into the halo of their host galaxies 



e.g.,|Marlowe et al. 1995; Martin 1998; Bomans 200l| [van" 



Eymeren et al. |20Q7 ) . This leads to the question of what the 
fate of the gas is. One scenario is that the gas cools down 
in the halo and eventually falls back onto the galactic disc 
(outflow, galactic fountain scenario, |Shapiro fc Field|1976| ). 
However, it might also be possible that the gas escapes from 
the gravitational potential by becoming a freely flowing wind 
(galactic wind). The detection of large amounts of hot gas 
in the intergalactic medium (IGM) and the generally low 
metal content of dwarf galaxies support this scenario. 

Hydrodynamic simulations modelling superbubble 
blow-out in dwarf galaxies ( |Mac Low &; Ferrara| [l999) pre- 
dict that at least a part of the hot gas has enough kinetic 
energy to leave the gravitational potential of its host galaxy 
and to enrich the IGM. The relatively low escape velocities 
of dwarf galaxies should facilitate the removal of substan- 
tial amounts of gas ( |Larson|[l974| ). According to |Ferrara fe| 



Tolstoy (2000) galactic winds occur in dwarf galaxies with 
gas masses up to 10 9 M . However, the fate of the gas does 
not only depend on the mass of the host galaxy and there- 
fore its gravitational potential, but it also strongly depends 
on the morphology of the ISM distribution. For a spherical 
distribution, the ISM seems to be more resistant to ejection 
( Silich fe Tenorio-Tagle(200l ). 

A recent detailed kinematic study of the neutral and 
ionised gas in the two irregular dwarf galaxies NGC 2366 



and NGC 4861 revealed several outflows in each galaxy (van 
Eymeren et al. 2009a b). The measured expansion velocities 
were of the order of 20 to 50kms _1 , which is, in compari- 
son with the escape velocities of the host galaxies, too low 
to allow the gas to leave the gravitational potential. This 
result confirms earlier observations: no convincing case of a 
galactic wind in local dwarf galaxies has been reported so 
far (Bomans 2005). Note that galaxies like M82 are not typ- 
ical dwarf galaxies as they show strong starbursts, are more 
luminous and have a higher mass. 

Altogether, this shows that despite all the detailed ob- 
servational studies and simulations, we do not fully under- 
stand what is going on in these galaxies. Galactic winds 
are still thought to be necessary ingredients for their forma- 
tion and evolution, but direct evidence seems to be difficult 
to get. In order to improve our knowledge about the pro- 
cesses happening in dwarf galaxies, we performed a multi- 
wavelength study of altogether four irregular dwarf galax- 
van Eymeren] |2008) . The results for NGC 2366 and 



NGC 4861 have already been published (van Eym eren et al 
2009a b, see also above). We here concentrate on the two 
remaining galaxies IC4662 and NGC 5408. We obtained Hi 



line data as well as optical images and spectroscopic data. 
Some basic properties of both galaxies are given in Table [I] 
IC4662 (HIPASS J1747-64) is classified as a barred 
irregular galaxy of Magellanic type (IBm). Its distance of 
Dtrgb = 2.44 Mpc was obtained by Karachents ev et al.| 
(2006) and hence makes it the nearest known representa- 
tive of blue compact dwarfs ( |Karachentsev et al.||2006| ). It 
seems to be a rather isolated galaxy, belonging to no known 
de Vaucouleurs Jl975l describes IC4662 as 



groups. 



fore- 



ground galaxy in the direction of the NGC 6300 group. 

NGC 5408 (HIPASS J1403-41) is classified as an 
IB(s)m galaxy. It was first studied by Bohuski et al. (1972) 
who found that its nucleus consists of several bright Hn 
regions and appears to be undergoing a violent burst of 
star formation. As the galaxy reveals an ultra-luminous 
X-ray source, NGC 5408 X-l, very close to the main Hn 
regions, it was a popular object to study over the last 
decade (e.g., Soria et al. 2004 2006 Strohmayer et al. 



2007 Lang et al. 2007 Kaaret & Corbel 2009). This 



X-ray source has recently been argued to harbour an 



intermediate-mass black hole (Strohmayer & Mushotzky 



[2009]). The dista nce of NGC 5408 of £>trg b = 4.81 Mpc 
was obtained by Karachentsev et al. (2002). Its position 



on the sky puts NGC 5408 in the Centaurus A group. The 
closest known neighbour appears to be ESO 325-G7001 at 
a projected distance of 208'. 

This paper is organised as follows: in Sect. 2 we give 
an overview over the observations and the data reduction. 
In Sect. 3 we describe and compare the morphology of the 
neutral and ionised gas. Section 4 contains the kinematic 
analysis of both gas components including a search for ex- 
panding gas. The results are discussed in Sect. 5, which is 
followed by a summary in Sect. 6. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical data 

2.1.1 Imaging 

Deep i?-band and Ha images are available for both galaxies. 
Using the ESO Multi-Mode Instrument (EMMI) attached to 
the ESO New Technology Telescope (NTT), we obtained a 
600 s exposure of IC 4662 in i?-band and a 1800 s exposure in 
Ha. A 600 s i?-band image of NGC 5408 as well as a 1200 s 
Ha image, both observed with the ESO Faint Object Spec- 
trograph and Camera (EFOSC) attached to the ESO 3.6m 
telescope, have been taken from the ESO archive. The data 
reduction was performed using the software package IRAF 
(Tod y|1993| , and included standard procedures of overscan- 
and bias-subtraction as well as a natfield correction. Addi- 



tionally, we removed cosmic rays by running the IRAF ver- 
sion of L. A. Cosmic ( van Dokkum|200T ). In order to get the 
pure Ha line emission, we first scaled the flux of the stars 
in both the continuum and the Ha images and subsequently 
subtracted the continuum image from the Ha image. The 
seeing was l'/l and 0'/8 during the observations of IC4662 
and NGC 5408 respectively. We used adaptive filters based 
on the H-transform algorithm (Ric hter et al.|1 991) to stress 
weaker structures in the Ha images and to differentiate them 
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Table 1. The basic properties of IC 4662 and NGC 5408. 



Optical name 
HIPASS name 


IC 4662 
HIPASS J 1747-64 


NGC 5408 
HIPASS J1403-41 


Ref. 


opt. centre: 






(1) 


a (J2000) 


17 h 47 m 08.8 s 


14 h 03 m 20.9 s 




S (J2000) 


-64° 38' 30" 


-41° 22' 40" 






2.44 


4.81 




^opt L^iii D J 


336 ± 29 


537 ± 33 


(1) 


type 


IBm 


IB(s)m 


(1) 


nnt HiflTTiptpr F'l 

\ 1 1 y L . v 1 1 1 1 1 v 1 v 1 


3.0 x 1.6 


2.6 x 1.6 


(4) 


nnf nlnmpfpr rifTM"" 


2.1 x 1.1 


S 6 x 2 2 




inclination [°] 


58 


52 


(4) 


position angle [°] 


1 n/i 


DZ 


w 


A B [mag] 


0.303 


0.298 


(5) 


ra B [mag] 


12.33 ±0.09 


12.59 ±0.09 


(4) 


B-V 


0.41 


0.56 




M B [mag] 


-14.91 ±0.1 


-16.12 ±0.09 




L B [10 9 L ] 


0.14 ±0.01 


0.44 ±0.03 




-uhi [kms -1 ] 


302 ±3 


506 ±3 


(6) 


i; L g [kms -1 ] 


153 


314 


(6) 


w 50 [kms -1 ] 


86 


62 


(6) 


^20 [kms -1 ] 


133 


112 


(6) 


F H i [Jykms -1 ] 


130.0 ±12.0 


61.5 ±6.7 


(6) 


M m [10 8 M ] 


1.83 ±0.17 


3.36 ±0.37 





Note: The blue luminosity is calculated using a solar B-band magnitude of 5.48 mag. - References: (l) |de Vaucouleurs et al.|(1991|), (2) 
|Karachentsev et al.| i f2006|) , (3) |Karachentsev et al~] ( |2002| ) 7 (4) |Lauberts fe Valentjjn] ( |1989| ) [ESO Uppsala], (5) |Schlegel et al. \1998\ , 
(6) |Koribalski et al.|p004| [HIPASS BGC]. 



Table 2. Imaging - some observational parameters. 



Parameter [Unit] 



IC4662 



NGC 5408 



Telescope 
Instrument 
Filter 
Date 

Exp. Time [s] 
Spatial res. [ /; 
Airmass 
Seeing ["] 



ESO NTT 
EMMI 
608 (R), 596 (Ra) 
07.05.91 
1x600, 1x1800 
0.36 
1.23 
1.1 



ESO 3.6m 
EFOSC 
642 (R), 692 (Ha) 
10.07.02 
1x600, 1x1200 
0.16 
1.02 
0.8 



from the noise. All images are displayed in Fig. [T] Some ob- 
servational details are given in Table [2] 



2.1.2 Medium-resolution long-slit spectroscopy 

We used archival medium-resolution long-slit spectra of 
NGC 5408, which have a spectral resolution of about 
60 kms -1 , as measured from the night sky lines. As the 
target of these observations was the ultra-luminous X-ray 
source NGC 5408 X-l, all spectra were obtained at almost 
the same position. Therefore, we only analyse one slit posi- 
tion as a representative. The seeing was 0'!9. 

We performed further medium-resolution long-slit spec- 
troscopy not only of NGC 5408, but also of IC 4662 from May 



1st to May 3rd 2006. As we aimed for very deep spectra and 
were mainly interested in the Ha line, only the red part of 
EMMI was read out. We used grating # 7 centred on the Ha 
line with a wavelength coverage of 1300 A and a dispersion 
of 0.41 A/pix. With a slit width of l" we achieved a spectral 
resolution of 112 kms -1 , again as measured from the night 
sky lines. The pixel size is ; /332 and the seeing was about 
l" . Spectra of a thorium comparison lamp were taken for 
the wavelength calibration. 

We obtained three spectra each at different positions 
across the two galaxies. The slit positions were chosen to 
intersect prominent shell and filamentary structures as visi- 
ble on the Ha images. Some details of the observations are 
given m Table [3] 

The data reduction of all spectra (including the archival 
one) was performed by us using IRAF and included 
overscan- and bias-subtraction, flatfield correction and a 
wavelength calibration. A geometric distortion correction 
was not necessary because the deviations were smaller than 
one pixel. Again, cosmic rays were removed by running L. A. 
Cosmic. We also performed a background subtraction to re- 
move the contamination by the night sky lines. 

The spectra and their analysis are presented in Sect. |4.4) 
They were binned by three pixels in spatial direction, which 
corresponds to about l" in order to match the seeing (see 
above). At positions of very weak emission, ten pixels were 
summed up to improve the signal to noise ratio (S/N). 
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30 s 25 s 14 h 3 m 20 s 15 s 30 s 25 s 1 4 h 3 m 20 s 1 5 : 

RA (2000.0) R.A. (2000.0) 



Figure 1. Optical imaging of the Magellanic irregular dwarf galaxies IC4662 (top) and NGC 5408 (bottom). Left panel: i?-band image; 
right panel: continuum-subtracted Ha image. The images of IC 4662 were obtained with the ESO NTT, the images of NGC 5408 with 
the ESO 3.6m telescope. In order to locate the positions of the main star clusters and the main Hll regions respectively, the high intensity 
areas in both images are highlighted in white. 



Table 3. Long-slit spectroscopy - some observational details. 



Object 


Telescope / Instrument 


Date 


Slit No. 


Exp. Time 


PA 










[min] 


[°] 


NGC 5408 


ESO NTT / EMMI 


01.05.06 


la 


5x45 


54 






02.05.06 


2 


4x45 


6 






03.05.06 


3 


4x45, 1x21.5 


296 






17.05.93 


lb 


3x20 


80 


IC 4662 


ESO NTT / EMMI 


02.05.06 


1 


4x45 


330 






03.05.06 


2 


3x45 


56 






02.05.06 


3 


4x45 


70 



Note: Slit lb of NGC 5408 denotes the archival spectrum. The position angle PA is defined to increase counter-clockwise with north 
being 0°. 
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2.2 ATCA Hi synthesis observations 

H I line observations of IC 4662 and NGC 5408 were obtained 
with the Australia Telescope Compact Array (ATCA), and 
are part of the "Local Volume H I Survey" (LVHI^ Korib- 
alski 20081 



Table 4. ATCA Hi data - imaging properties. 



Koribalski et al.|2010 ). 



The first frequency band was centred on 1418 MHz with 
a bandwidth of 8 MHz, divided into 512 channels. This gives 
a channel width of 3.3 km s -1 and a velocity resolution of 
4kms _1 . The ATCA primary beam is 33'6 at 1418 MHz, 
which is sufficient to fully map both galaxies and their sur- 
roundings. 

In order to ensure excellent uv- coverage and sensitiv- 
ity to large-scale structures, IC4662 was observed for a full 
synthesis (12 h) in each the EW367, 750A, and 1.5D config- 
urations. The lower resolution data of NGC 5408 (375 and 
750D array) were taken from the archive and complemented 
by high-resolution data within LVHIS (1.5 A array). 

The data reduction was carried out with the MIRIAD 



software package (Sault et al. 1995) using standard proce- 
dures, including flux- and phase calibrations. Using a first 
order fit to the line- free channels in the Hi data set, the 
20 cm radio continuum was separated from the H I emission. 
We created two different sets of data cubes, a low-resolution 
cube to detect the extended emission and a high-resolution 
cube to resolve the inner structure and to make the H I data 
comparable to the optical data. For the low-resolution cubes 
of both galaxies, we excluded the longest baselines, which 
are all baselines to the distant antenna 6 (CA06). The low- 
resolution data of IC 4662 were made using 'natural' weight- 
ing of the uv-d&ta, in the velocity range covered by the Hi 
emission in steps of 4kms _1 . The low-resolution data of 
NGC 5408 were made using 'robust' weighting. The high- 
resolution data of both galaxies were created by including 
all baselines to CA06 and using the same weighting as for 
the low-resolution cubes. The resulting beam sizes are given 
in Table H 

Uybeam -1 corresponds to an Hi column density of 
2.55 x 10 20 atoms cm" 2 (IC4662, low-res.), 2.64 x 10 21 
atoms cm -2 (IC4662, high-res.), 2.89 x 10 20 atoms cm -2 
(NGC 5408, low-res.), and 2.77 x 10 21 atoms cm" 2 
(NGC 5408, high-res.). The moment maps (integrated in- 
tensity map, intensity-weighted mean velocity field, and the 
velocity dispersion) were created from the Hi data cubes 
by first isolating the regions of significant emission in ev- 
ery channel and afterwards clipping everything below a 2.5a 
threshold. This final step of the data reduction process and 
the subsequent analysis of the H I data were performed with 



The Groningen Image Processing System (GIPSY, van der 
Hulst et al.||1992| ), complemented by some IRAF tasks. 



1 LVHIS is a large project that aims to provide detailed H I inten- 
sity maps, velocity fields and 20 cm radio continuum observations 
for a complete sample of nearby, gas-rich galaxies belonging to the 
Local Volume (LV), a sphere of radius lOMpc centred on the Lo- 
cal Group. The ATCA observations include all LV galaxies which 
were detected in HIPASS and reside south of 5 « -30°. Further de- 
tails can be found at http://www.atnf.csiro.au/research/LVHIS/. 



Object 


IC 4662 


NGC 5408 


Low-resolution: 






Weighting 


natural, -CA06 


robust, -CA06 


H I synthesised beam 


70" x62" 


71" x 54" 


High-resolution: 






Weighting 


natural, +CA06 


robust, +CA06 


H I synthesised beam 


21" x 20" 


20" x 20" 



3 GENERAL MORPHOLOGY 

In order to get a first impression of how the stars and the 
gas are distributed, we compared the optical and Hi mor- 
phologies of IC4662 and NGC 5408. 



3.1 IC4662 

The i?-band image of IC 4662 (see Fig. [I] upper left panel) 
reveals a slightly elongated, box- like shape with a position 
angle of 104° (Table [I]). The Ha image (Fig. [I] upper right 
panel) shows a large H n region complex with an Ha mini- 
mum in the south-eastern part of the stellar light distribu- 
tion, and extended diffuse, filamentary gas structures out- 
side the stellar disc to the north-east. The central two star 
clusters are offset by about 350 pc to the north-east from the 
optical centre (see Table [I]), which probably explains the ex- 
tended Ha emission in this area. They coincide with the two 
brightest H II regions visible on the Ha image. Additionally, 
a star cluster with associated Ha emission that appears to 
be detached from the main body of the galaxy can be seen. It 
is located 1'5 or 1.1 kpc to the south-east of the centre. Our 
deep Ha image shows an Hn region at the same position 
that is connected to the main complex by small, compact 
Hn regions and diffuse filamentary gas structures. We will 
come back to this feature in Sect. [5751 

The channel maps of the low-resolution Hi data (see 
Fig. [2]) show that the emission is distributed over a velocity 
range from roughly 200 to 400kms _1 , i.e., over 200kms _1 . 
The neutral gas seems to consist of two systems, an inner 
system with a position angle of about 135° and an outer part 
with a position angle of roughly 45°, i.e., perpendicular to 
the inner system. The integrated Hi intensity map (upper 
left panel of Fig. [3]) confirms this feature: it shows an inner 
high column density system that is perpendicular to the 
outer low column density system. As displayed on the lower 
right panel of Fig. [3] the inner system coincides with the 
optical extent of IC 4662 including the southern H 11 region. 
We can see an additional Hi maximum to the north-west 
of the galaxy, which has no optical counterpart. The Hi 
outer diameter is given in Table [5] Note that the neutral 
gas is extended by a factor of about six in comparison to 
the ionised gas, which is among the largest ratios of H 1 over 
Haever detected. We measured the flux density to be Fm — 
123Jykms _1 , from which we derived a total Hi mass of 
Mm = 1.7 x 10 8 M (see Table |HJ. This is within the errors 
of the HIPASS value of F m = (130.0 ± 12.0) Jy kms" 1 (see 
Table [TJ, which means that our interferometric observations 
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Table 5. ATCA Hi properties. 



Optical name 


IC 4662 


NGC 5408 


HIPASS name 


J1747-64 


J1403-41 


dynamic centre: 






oi(J2000) 




14 h 03 m 21.9 s 


5(J2000) 




-41° 22' 03" 


<u sys [kms -1 ] 


302 


502 






62 


PA [°] 




302 


F H i [Jykms- 1 ] 


123 


63 


M m [10 8 M ] 


1.7 


3.4 


H I diameter ['] 


15.0 x 11.5 


10.0 x 6.3 


H I diameter [kpc] 


10.6 x 8.2 


14.0 x 8.8 


Hi / opt. ratio 


6 


4 


v TO t [kms -1 ] 




53 


< a > [kms -1 ] 


11/20 


17 


(Jpeak [kms -1 ] 


37 


27 


r*HI,max [kpc] 


6 


9.5 


M dyn [10 9 M ] 




6.2 



Note: As the velocity field of IC 4662 is very distorted, we did 
not perform a tilted-ring analysis and therefore could not 
calculate any of the kinematic parameters. The systemic 
velocity was estimated from the global Hi profile of IC4662. 
The lower value of < a > refers to the outer system of IC 4662, 
the higher value to the inner, perpendicular system. 

recover most of the emission detected on single-dish data of 
IC4662. 

We also show the Hi intensity distribution at higher 
spatial resolution (21" x 20") for a better comparison with 
the optical data. In the upper left panel of Fig. [4] the high- 
resolution Hi intensity contours are overlaid in white onto 
the continuum-subtracted Ha image. Additionally, the in- 
nermost contours of the low-resolution H I intensity distribu- 
tion are displayed in black. As the high-resolution contours 
show, the central elongated maximum is split into several 
smaller, point-like maxima. Interestingly, all these maxima 
are clearly offset from the most prominent Hn regions. 

3.2 NGC 5408 

Like IC4662, NGC 5408 reveals a slightly elongated stellar 
disc with a position angle of 62°, which leads to a box-like 
appearance in the optical (see Fig. [I] lower left panel). The 
main star-forming knots are located along a chain running 
from the centre to the main star formation area in the south- 
west. In Ha, the galaxy is dominated by a bright complex of 
Hll regions in the south-west (Fig.[l] lower right panel) that 
coincides with the brightest star clusters. In the eastern part, 
the ionised gas becomes more patchy with a few holes in the 
distribution and several filaments emanating from the disc. 
Indeed, it seems that the eastern part of the galaxy ends in a 
superbubble as suggested by the shell-like structures visible 
in the north- and south-east and the big hole in between. 
Diffuse gas can also be seen across the southern part of the 
galaxy. 

Figure [5] shows the low-resolution H I channel maps of 
NGC 5408. We detected emission over 130 km s -1 , from 430 
to 560kms _1 . In the inner parts the emission is elongated 



along the east-west axis (position angle of 90°), whereas in 
the outer parts it is more extended to the north-west and 
south-east (position angle of 135°). The integrated Hi in- 
tensity distribution looks fairly symmetric and reveals two 
point-like maxima in the centre (see Fig. [5] and upper left 
panel of Fig. [6}. A comparison with the ionised gas distri- 
bution shows that the western Hi intensity maximum co- 
incides with the centre of the optical emission in NGC 5408 
(see Fig. [6] lower right panel). On larger scales, the distribu- 
tion of the neutral gas seems to be warped in the outer parts, 
which explains the two orientations seen in the channel maps 
(Fig. [5}. In comparison to the ionised gas, the neutral gas 
component is four times more extended, which is again an 
unusually high ratio (for the H 1 diameter see Table 5 . We 
measured the flux density to be Fhi = 63 Jy km s~ , from 
which we derived a total Hi mass of Mm = 3.4 x 10 8 Mq 
(Table [5). This is again within the errors of the HIPASS 
value of Fhi = (61.5 =b 6.7) Jykms -1 . 

Again, we overlaid the high-resolution H 1 intensity con- 
tours (20" x 20") in white over the continuum-subtracted 
Ha image (Fig. [4] upper right panel). The distribution shows 
two prominent maxima which coincide with the point-like 
maxima detected on the low-resolution H I map (black con- 
tours). The western extended maximum coincides with the 
main Hll region complex in NGC 5408, although it is not 
fully centred on the brightest Hll region. It is connected 
via a bridge with the eastern maximum which is clearly off- 
set from the southern shell-like structure mentioned above. 
Additionally, the high-resolution Hi intensity distribution 
reveals a third, weaker maximum in the north, close to, but 
offset from the northern shell-like structure. In agreement 
with the distribution of the ionised gas, the H 1 column den- 
sity between the shell-like structures is very low. This con- 
firms our findings in the optical data, i.e., the existence of a 
superbubble in the eastern part of NGC 5408. 



4 KINEMATIC ANALYSIS 

The comparison of the optical and H 1 morphology in IC 4662 
and NGC 5408 showed that on low-resolution scales, the op- 
tical extent coincides very well with the central Hi maxi- 
mum. However, the high-resolution H 1 maps reveal that the 
peaks in the neutral gas distribution are typically offset from 
the Hn regions. On larger scales, the optical appearance of 
both galaxies is defined by an alignment which is perpen- 
dicular to the extended H 1 distribution (IC 4662) or differs 
by at least 45° (NGC 5408). 

In this section, we analyse the kinematics of the two 
galaxies in order to study the properties of both their neu- 
tral and ionised gas components. All given velocities are he- 
liocentric velocities measured along the line of sight. 

4.1 The distorted H 1 kinematics of IC 4662 

The Hi velocity field of IC4662 (Fig. [3] upper right panel; 
see also Fig. [4] middle row, left panel) is very disturbed. 
The overall velocity gradient runs from the north-east with 
velocities of 220kms _1 to the south-west with velocities of 
380kms _1 . The direction of the rotation of the neutral gas 
seems to change from the western part of the galaxy to the 
eastern part, causing a sudden change of the position angle 
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Figure 2. Low-resolution Hi channel maps of IC4662 (contours) superposed on our i?-band image. For display purposes we show the 
maps with a velocity resolution of 12 kms -1 . The original channel spacing is 4 kms -1 . The noise is about 0.9 mjy beam -1 . Contours are 
drawn at -2.5 (-3cr), 2.5 (3a), 5, 10, 20, 40, 80, and 160 mJy beam -1 . The synthesised beam of 70 // x62 // (natural weighting, excluding 
CA06) is displayed in the lower left corner of the first channel map. The optical centre of the galaxy is marked by a white cross in the 
same channel map. The corresponding heliocentric velocities are indicated in the upper right corner of each channel. 2' correspond to 
about 1.4 kpc. 



of almost 90°. This change is possibly connected with the 
two perpendicular systems as discussed in Sect. |3.1| The 
H I velocity dispersion map (lower left panel of Fig.[3] over- 
laid are the H I intensity contours; see also Fig. [4] lower 
left panel) shows a dispersion of 20 kms -1 in the inner and 
11 kms -1 in the outer parts, except for two maxima with al- 
most 40 kms -1 east of the extended Hi intensity maximum. 

The large-scale twist in the velocity field makes it im- 
possible to derive a rotation curve for IC 4662. We extracted 
spectra from the Hi data cube to get a better understand- 
ing of the kinematics of the neutral gas. They are shown 



in Fig. |7| each box represents the sum of all spectra within 
70" x 70" (corresponding to one beam size). Stars denote 
spectra with a flux scale which is four times larger than 
that of the other spectra. It is obvious that across the whole 
distorted area, the Hi line is often split into two, some- 
times clearly separated components. The difference of about 
70 kms -1 is always the same, although the absolute veloci- 
ties vary according to the rotation of IC4662. 

As already mentioned in Sect. |2.2| we created the veloc- 
ity field using the intensity- weighted mean (iwm). However, 
as soon as the velocity distribution of the H I emission is not 
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Figure 3. Hi moment maps of IC4662, combining three arrays, but excluding baselines to CA06, which leads to a synthesised beam 
of 70 ;/ x 62 /; (natural weighting). Top left: the Hi intensity distribution. Contours are drawn at 0.06 (3<x), 0.16, 0.4, 0.8, 1.6, 3.2, and 
6.4 Jy beam -1 kms -1 where Uybeam -1 corresponds to a column density of 2.55 x 10 20 atoms cm -2 . Top right: the Hi intensity- 
weighted mean velocity field. Contour levels range from 220 to 380 kms -1 in steps of 10 kms -1 . The systemic velocity of 302 kms -1 , 
estimated from the global Hi profile, is marked in bold. Bottom left: the Hi velocity dispersion. Overlaid are the same Hi intensity 
contours as mentioned above. Bottom right: continuum-subtracted Ha image. The spatial resolution was smoothed to fit the resolution 
of the Hi data. Overlaid in black are again the Hi intensity contours. Overlaid in grey are the Hi velocity dispersion contours at 
25 kms -1 . 



symmetric, the velocity values are biased to the longest tail 
of the velocity profiles ( de Blok et al.||2008 ). Therefore, one 
might argue that the iwm velocity field does not represent 
the true velocities. A recently more popular approach is to 
fit the line profiles with Gauss-Hermite polynomials, which 
allows to define the peak velocities more accurately. This be- 
comes important when components overlap. In the IC 4662 
spectra, the two components are clearly separated with one 
of them dominating the H I profile. Hence, the difference be- 
tween the iwm and the Hermite velocity field should be neg- 
ligible. Nevertheless, we created the Hermite velocity field 
by fitting Gauss-Hermite h3 polynomials to all line profiles 



using the GIPSY task xgaufit (Fig. Al middle panel). The 
right panel of Fig. |A1| shows the residuals after subtracting 
the Hermite velocity field from the iwm. As we expected, the 
residuals are mostly close to zero, except for an area in the 
north-east. A look at the spectra (Fig. |7| reveals that the 
S/N in this area is quite low, which can often not be han- 
dled properly by the GIPSY moments task which we used to 
create the iwm velocity field. As the large-scale distortion is 
still very pronounced in the Hermite velocity field, we think 
that it is not a consequence of using the 'wrong' method 



to create the velocity field, but that it is a real feature in 
IC4662. 



4.2 The H I velocity field of NGC 5408 

NGC 5408 shows a symmetric H I velocity field (Fig. [6] up- 
per right panel; see also Fig. [4] middle row, right panel) 
with a velocity gradient running from the south-east with 
velocities of 450 kms -1 to the north-west with velocities of 
550 km s _1 . The gradient is smooth, even at both ends of the 
disc where the Hi distribution looks warped. It even seems 
that it better follows the large scale structure than the in- 
ner east- west alignment. A local, quite prominent distortion 
of the velocity field can be seen south of the Hi intensity 
maxima. The H I intensity distribution, however, appears to 
be smooth and undisturbed (see upper left panel of Fig.|6|. 
The H I velocity dispersion (lower left panel of Fig. [6] over- 
laid are the Hi intensity contours; see also Fig.^] lower right 
panel) shows two peaks at 25 kms -1 , which are located on 
either side of the spur in the velocity field. 

Again, we extracted the spectra from the H I data cube 
in order to analyse the profiles (see Fig.j8|. Each box equals 
a beam size of 60" x 60". Stars denote spectra with a flux 
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Figure 4. High-resolution Hi moment maps of IC4662 (left) and NGC 5408 (right). Upper row: a comparison of the Ha and Hi line 
emission in the inner discs of both galaxies. The continuum-subtracted Ha emission is shown as grey-scale, while the black and white 
contours indicate the H I emission at low and high angular resolution respectively. The high-resolution contours of IC 4662 are drawn 
at 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, and 1.6 Jy beam -1 kms -1 where one beam corresponds to 21" x 20"; the ones of NGC 5408 at 1.2, 
1.6, 2.0, 2.4, 2.8, 3.2, 3.6, and 4 Jybeam -1 kms -1 where one beam corresponds to 20" x 20". Middle row: Hi velocity fields. Contour 
levels are given in steps of 10 kms -1 and range from 260 to 340 kms -1 (IC4662) and from 480 to 530 kms -1 (NGC 5408). The systemic 
velocities (302 kms -1 and 502 kms -1 respectively) are marked in bold. Lower row: Hi velocity dispersion maps. The high-resolution 
H I intensity contours are overlaid as in the top panels. 
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Figure 5. Low-resolution Hi channel maps of NGC 5408 (contours) superposed on the i?-band image. For display purposes we show the 
maps with a velocity resolution of 12kms — 1 . The original channel spacing is 4kms _1 . The noise is about lmjybeam - 1 . Contours are 
drawn at -3 (-3<r), 3 (3a), 5, 10, 20, 40, 80, and 160 mJy beam -1 . Otherwise the same as in Fig. [3] The synthesised beam is 71" X 54" 
(robust weighting, excluding CA06). 2' correspond to 2.8 kpc. 



scale which is four times larger than that of the other spec- 
tra. In general, the Hi profiles look quite simple. The lines 
are broad with FWHM of 35 to 50kms _1 , but often also 
symmetric with no second component clearly standing out. 
The high velocity dispersion in the south coincides with an 
area of very low S/N, which is probably again due to the in- 
ability of the moments task to properly handle noisy spectra. 
However, at the eastern edge we detected two components 
which are separated by about 40kms -1 . 

We tested the quality of the iwm velocity field by sub- 
tracting the Hermite velocity field: the residuals are very 
small except for the area of high velocity dispersion. As in 



the case of IC4662, this can be explained by the low S/N 
and the line-splitting. For a direct comparison we refer the 



reader to van Eymeren et al. ( 2009 ) where the H I kinemat- 



ics of several dwarf galaxies (among them NGC 5408) have 
been analysed by using the Hermite velocity fields. 



4.3 The H I rotation curve of NGC 5408 

In order to make predictions about the fate of the gas, we 
need to know some of the kinematic parameters like the in- 
clination and rotation velocity of NGC 5408. Therefore, we 
derived a rotation curve by fitting a tilted-ring model to the 
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Figure 6. Hi moment maps of NGC 5408, combining three arrays, but excluding all baselines to CA06, which leads to a synthesised 
beam of 71 /x x 54 /; (robust weighting). Top left: the Hi intensity distribution. Contours are drawn at 0.2 (3<x), 0.4, 0.8, 1.2, 2.0, 3.2, 
and 4.0 Jy beam -1 kms -1 where Uybeam - 1 correspond to a column density of 2.89 x 10 20 atoms cm -2 . Top right: the Hi intensity- 
weighted mean velocity field. Contour levels range from 450 to 550 kms -1 in steps of 10 kms -1 . The systemic velocity of 502 kms -1 , 
as calculated from fitting tilted-rings to the velocity field, is marked in bold. Bottom left: the Hi velocity dispersion, overlaid are the 
same Hi intensity contours as mentioned above. Bottom right: continuum-subtracted Ha image. The spatial resolution was smoothed 
to fit the resolution of the Hi data. Overlaid in black are again the Hi intensity contours. Overlaid in grey are the Hi velocity dispersion 
contours at 20 kms -1 . 



observed velocity field (GIPSY task rotcur, Begeman 1989). 
Initial estimates for the relevant parameters (systemic veloc- 
ity, the coordinates of the dynamic centre, the inclination, 
and the position angle) were obtained by interactively fit- 
ting ellipses to the H I intensity distribution with the GIPSY 
task ellfit. The width of the rings was chosen to be half the 
spatial resolution, i.e., 31". In order to get the most precise 
values, three different approaches were made by always com- 
bining receding and approaching side. The upper left panel 
of Fig. [9] shows the resulting curves: first, we kept the ini- 
tial estimates fixed in order to derive the rotation velocities 
(green (light grey) triangles). The best-fitting values for all 
parameters were derived by allowing only one parameter to 
vary with radius while keeping the remaining ones fixed. As 
no significant variation of any of the parameters over radius 
was noticed, we derived an average value (given in Table [5}. 
The black triangles represent the resulting rotation veloci- 
ties. The error bars indicate the values of the receding (top) 
and approaching side (bottom) when treated separately. In 



order to reproduce the result of the iterative approach, the 
so-derived parameters were all left free (red (dark grey) tri- 
angles). As the filling factor of the rings drops significantly 
beyond a radius of 320", we did not fit the outer parts. That 
means we also miss most of the distortion at the southern 
edge of the galaxy. 

NGC 5408 has a slowly rising rotation curve coming to 
a plateau at a radius of 250". Receding and approaching side 
show a very similar kinematic behaviour up to a radius of 
320", confirming the impression of an evenly rotating galaxy. 
The initial estimates already defined the H 1 kinematics quite 
well as the green and black triangles are in very good agree- 
ment. The 'left-free' approach is also in good agreement with 
the other two approaches, except for some scatter between 
0" and 100". This is still far away from the disturbed region, 
which can therefore not explain the scatter. 

The iterative approach results in an inclination angle 
of 62°, a position angle of 302°, and a systemic velocity 
of 502 kms -1 (see also Table pk. The systemic velocity is 
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Figure 7. H I velocity profiles of IC 4662. We show the low-resolution H I velocity field (see Fig. [5] upper right panel) with the H I profiles 
overlaid. Each box represents the sum of all spectra within roughly one beam size (70" x 70 ;/ ). Stars denote spectra with a flux scale 
which is four times larger than that of the other spectra. 



slightly lower than the HIPASS velocity of 506 km s" 1 (Ko 



ribalski et al. 2004). In comparison to the parameters derived 



from optical data, the inclination of the neutral gas is higher 
by 10°. The value of the position angle of the neutral gas is 
dominated by the extended emission. As the optical disc is 



differently aligned (see Sect. 3.2), the Hi position angle of 
302° differs significantly from the optical position angle of 
60°. 



In order to prove the reliability of the calculated values, 
a model velocity field was created by using the best-fitting 
parameters (Fig. [9] lower left panel) and afterwards sub- 
tracted from the original velocity field (upper right panel). 
The residual map is shown on the lower right panel of Fig. [9] 
with the 3cr contour from the low-resolution Hi intensity 
map and the systemic velocity contour overplotted. It is ob- 
vious that the original velocity field is very well described by 
the derived kinematic parameters. The residuals are gener- 
ally at zb 5kms -1 , except for one region in the south, which 
shows residuals of about 15kms _1 . This region adjoins the 
already mentioned distortion in the velocity field and the 
high velocity dispersion. We examine this issue in more de- 
tail in Sect. l4A2l 



4.4 A detailed kinematic analysis of the H I and 
Ha line data 

We now want to study the H I and Ha line profiles of both 
galaxies in more detail and to compare the peak velocities 
of the neutral and ionised gas. Therefore, we performed a 
Gaussian decomposition by interactively fitting the H I and 
Ha line emission (IRAF task splot). Only detections above 
a 3cr threshold were considered. 

Figures [7] and [8] have already shown that the Hi profiles 
of both galaxies are sometimes split into at least two com- 
ponents. For a better comparison with the Ha images, we 



used the high-resolution H I data cubes (see Sect. 2.2 ) to per- 
form the Gaussian decomposition. As the extracted spectra 
are often quite noisy, especially in the outer parts of both 
galaxies, only the central parts are fitted. Figures [10] and [11] 
show the resulting maps: a map with the blue-shifted gas 
on the left, a map with the main component (component of 
highest intensity) in the middle and a map with red-shifted 
gas on the right. The velocities are averaged over 20" x 20" 
in case of IC 4662 and 18" x 18" in case of NGC 5408, which 
roughly corresponds to one beam size. For better visuali- 
sation, we overlaid the 3a Hi intensity contour from the 
low-resolution data and the 3a Ha intensity contour from 
the Ha images (Fig.[T]). 

We then performed a Gaussian decomposition of the Ha 
line profiles as measured from the long-slit spectra. We mea- 
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Figure 8. Hi velocity profiles of NGC 5408. We show the low-resolution Hi velocity field (see Fig. [6] upper right panel) with the Hi 
profiles overlaid. Each box represents the sum of all spectra within roughly one beam size (60" X 60" ). Stars denote spectra with a flux 
scale which is four times larger than that of the other spectra. 



sured the peak velocities in order to create position- velocity 



line is a superposition of at least two components (see, e.g., 



(pv) diagrams (Figs. [12] lower row and 13 middle left panel 
and lower row). The Ha velocities are indicated by + sym- 
bols. The corresponding Hi velocities (solid lines) were ex- 
tracted from the peak velocity fields of the main component 
(Figs. 10 and |11| middle panels). Both the Hi and the Ha 
velocities are not corrected for inclination. Figures [12] and 
[T3] also show the slit positions, plotted over the Ha image 
of each galaxy. Position in the pv diagrams is marked by a 
small circle; the arrows indicate increasing distance from 
in a positive sense. Examples illustrating the quality of the 
spectra and the analysis are given for both data sets in the 
upper right panel of Fig. [12] and in the upper and middle 
right panels of Fig. [13] 



At first glance, the Ha emission line is quite broad along 
all slits. Observations of the ionised gas in other dwarf galax- 
ies revealed out flows with expansion velocities of about 20 to 
SOkms -1 



(e.g. 



Schwartz &; Martin|2004[[van Eymeren et al.| 



2009a b). For the data used in our analysis, this implies that 



it will be difficult to separate blended lines because of the 
relatively low spectral resolution (60kms -1 for the archival 
spectrum of NGC 5408, 112kms _1 for all other spectra). 
Therefore, we fitted the profiles with only one Gaussian un- 
less a fit with two Gaussians gave reasonable values, i.e., 
FWHM h s > FW HMinstr- In many cases, a single Gaus- 
sian fit led to large residuals, which indicates that the Ha 



Fig. 12 upper right panel). 

In Fig. [14] we present the FWHMs, corrected for 
instrumental broadening FWFLMl orr = FWHM^ bs - 
FWHM? nstr , of all fitted lines vs. the integrated intensity. 
Different shades of grey indicate FWHMs measured in spec- 
tra from different slit positions. The black stars denote the 
data of the archival spectrum of NGC 5408. 

In the following subsections, we describe the kinematics 
of the neutral and ionised gas, point out peculiarities, and 
compare the kinematic behaviour of both gas components. 



4.4.1 IC 4662 

The high-resolution H I velocity maps of IC 4662 show that 
the neutral gas often consists of a single component that 



follows the overall rotation (Fig. 10). Red-shifted gas can 
be seen in the filamentary north-eastern part of the opti- 
cal disc of IC4662, which is offset from the main compo- 
nent by about 55kms _1 . We did not find any more outflow- 
ing neutral gas associated with the optical disc. However, 
we detected blue-shifted gas north-west and south-east of 
the optical disc with expansion velocities of about 40 to 
50kms _1 with respect to the main component, whose mor- 
phology gives the impression of some kind of bipolar outflow. 
We then used the Hi velocity values of the main com- 
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Figure 9. The Hi rotation curve of NGC 5408 created by fitting a tilted-ring model to the observed velocity field. Top left: different 
approaches for deriving the rotation curve. The black triangles represent the iterative approach, the error bars indicate receding and 
approaching sides. The green (light grey) curve was derived by taking the initial estimates and keeping them fixed, the red (dark grey) 
curve by taking the best-fitting parameters and leaving them free. Top right: the observed Hi velocity field (see Fig. [6}. Bottom left: 
the model velocity field created by taking the best-fitting parameters. Bottom right: the residual map after subtracting the model from 
the original velocity field. The 3a contour from the low-resolution H I intensity map and the systemic velocity are overplotted in black. 
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Figure 10. Gaussian decomposition of the H I line profiles of IC 4662, which were extracted from the high-resolution cube. We show the 
blue-shifted (left panel), main (middle panel) and red-shifted (right panel) components of the Hi velocities. The 3a Hi intensity contour 
from the low-resolution cube and the 3a Ha intensity contour from the Ha image (Fig. Q are overlaid in black for better visualisation. 
We averaged the velocities over 20" x 20", which roughly corresponds to one beam size. 
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Figure 11. Gaussian decomposition of the Hi line profiles of NGC 5408. The same as in Fig. ^0] We averaged the velocities over 
18 /; x 18 /; , which roughly corresponds to one beam size. 
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Figure 12. Gaussian decomposition of the Ha line profiles of IC 4662. Upper left panel: the continuum-subtracted Ha image. White 
arrows indicate the slit positions, white circles mark position 0. Upper right panel: an example spectrum (slit 3 at r = 110") which 
shows the flux in arbitrary units vs. wavelength. The black solid line indicates the original spectrum, the light grey line marks the Gaussian 
fit. The residuals after subtracting the Gaussian from the original profile are shown in dark grey. Note that the spectral resolution is about 
2.5A, which corresponds to 112 km s -1 . Lower panels: the pv diagrams of the Ha emission. The + symbols represent the Ha velocities 
obtained from the spectra, the solid line the H I velocities obtained from the high-resolution velocity map of the main component (see 
Fig. [To] middle panel). 



ponent as a reference value for the velocities of the ionised 
gas. We obtained optical spectra of IC4662 along three slit 
positions (Fig. 



12 



lower row). Slit 1 intersects the southern 
Hll region, the main Hn region complex, and the diffuse 
gas in between. South of the main Hn region complex, i.e., 
from — 160" to — 80", the ionised gas roughly follows the 
Hi velocities. In most cases, however, the subtraction of the 



fitted Gaussian from the observed profile reveals two resid- 
uals which are symmetrically distributed. The FWHM is 
measured to be 70kms -1 , which suggests a blue- and a red- 
shifted component with a velocity offset of 35kms -1 with 
respect to the Hi line. Between — 80" and —15", the ionised 
gas appears to be slightly red-shifted by 20 to 30 km s _1 with 
respect to the neutral gas. 
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Figure 13. Gaussian decomposition of the Ha line profiles of NGC 5408. Upper left panel: the continuum-subtracted Ha image. 
Again, all slit positions are indicated by white arrows. Position is marked by white circles. Upper and middle right panels: two 
example spectra (slit la at r = — 80", slit lb at r = ;/ ). Black solid lines indicate the original spectra, light grey lines mark the Gaussian 
fit in case of slit la and the sum of two Gaussian fits (dark grey lines) in case of slit lb). The residuals are shown in dark grey. Middle 
left panel: the pv diagram showing the Ha velocities extracted from the archival spectrum; lower panels: the pv diagrams of the three 
spectra that were obtained by us. The solid lines represent the Hi velocities extracted from the velocity map of the main component (see 
Fig. Ill] middle panel). 



Slit 2 intersects the northern part of the main H n region 
complex. The Ha emission detected along this slit follows 
the H I velocities quite well apart from a deviation of up to 
20kms _1 red-shifted between 50" and 80". The FWHM of 
the lines in this area is quite low in comparison to the other 
positions along the slit, and the residuals after subtracting 
a single Gaussian are negligible. That means that here, we 
directly detected expanding gas located in the north-east of 
the optical disc. 



Slit 3 intersects the central part of the main H II region 
complex. The Ha velocities observed along this slit are in 
good agreement with the Hi velocities from 50" to 150". We 
measured a FWHM of about 70kms -1 and detected signif- 
icant residuals after subtracting a single Gaussian (see, e.g., 
Fig. 12 upper right panel). Therefore, the ionised gas prob- 
ably consists of two components, a blue- and a red-shifted 
one, expanding with 35kms _1 with respect to the neutral 
gas. In the eastern filamentary area, i.e., from 0" to 50", the 
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Figure 14. The Gaussian FWHM of the Ha line, as measured 
from the long-slit spectra of both galaxies after correcting for 
instrumental broadening, vs. the integrated intensity. The small 
panels in the upper right corner of both diagrams show an en- 
largement of the FWHMs of the low intensity lines. 



pv-diagram shows a significant deviation of the Ha velocities 
from the H I velocities. While the Ha velocities still seem to 
follow the velocity gradient of the neutral gas between 50" 
and 150", the Hi velocities drop by about 40kms _1 . Look- 
ing at the velocity map of the main component of the neutral 



gas (see Fig. 10 middle panel), it becomes clear that this 
drop is part of the large-scale velocity field and coincides 
with the twist in the central parts where the velocity gradi- 



ent changes its direction (see also Sect. 3.1). That indicates 
that the neutral and ionised gas are kinematically decou- 
pled. Still, if we extend the gradient between 50" and 150" 
down to 0", the ionised gas seems to be red-shifted by 20 to 
30kms -1 . This red-shifted gas component was also detected 
in Hi (see above). 

Summarised, it can be said that the neutral and ionised 
gas are mostly in good agreement. The Ha residuals often 
imply the existence of two components, symmetrically dis- 
tributed around the H I reference value. There are a few ob- 
vious discrepancies, which might indicate outflowing ionised 



very disturbed Hi velocity field (see Sect. 5.1) 



44.2 NGC 5408 

Similar to IC 4662, the H I line profiles of NGC 5408 often re- 
veal a single component that follows the overall rotation. Co- 
inciding with the northern part of the optical galaxy, the H I 
shows an additional component, which is blue-shifte d w ith 
respect to the main line by about 30kms _1 (see Fig. 11). 

Again, we used the H I peak velocities of the main com- 
ponent as a reference value for the Ha velocities. We anal- 
ysed the kinematics of the ionised gas within NGC 5408 us- 
ing spectra provided by two different observing runs (see 



Sect. 2.1.2). The panels located in the lower row of Fig. 13 



show the pv diagrams of our own data, while the middle left 
panel represents the archival data. 

In comparison to IC4662, the velocities of NGC 5408 
seem to deviate more and stronger from the Hi reference 
value. Although the Ha line profile is generally well fitted by 
a single Gaussian, we detected more outflows as the outflow- 
ing gas often seems to dominate the intensity distribution. 

The centre of the Ha emission in the south-west har- 
bours a strong blue-shifted component, which has a velocity 
offset of about 60kms _1 in comparison to the neutral gas 
(see slit la). In slit lb (the archival data), there is a hint 
of even higher blue-shifted ionised gas (up to lOOkms -1 ). 
However, as only three measurements of this fast expanding 
outflow are available, we will not further consider it. Addi- 
tionally, slit lb shows a red-shifted component with a similar 
offset as seen in slit la of about 50kms _1 . This component 
is only visible in the higher resolution data, although a few 
measurements have been made in slit la as well (at —60"). 

The southern part of the shell structure in the east 
Sect ~ 



3.2) can kinematically be detected in slit lb, 2, 



(see 
and 

3. Slit lb, that covers the edge of the shell, shows a blue- 
shifted component with a velocity offset of about 20kms _1 
(at —80") with respect to the Hi reference value. In slit 2, 
which fully intersects the potential superbubble, the blue- 
shifted gas is offset by 30kms _1 at 30". Additionally, we see 
red-shifted gas at 20 km s _1 along the whole slit and a hint of 
an even higher red-shifted outflow of 70 km s _1 at 40". Slit 3, 
that also fully intersects the potential superbubble, but at 
a different angle, shows blue-shifted gas at 30", offset by 
50 km s _1 , and at 80", offset by 30 km s _1 . The latter outflow 
coincides with the northern part of the shell. Red-shifted gas 
at 20 km s _1 is visible between the two blue-shifted outflows. 

As mentioned above, the neutral gas shows a hint of a 
blue-shifted component in the northern part of the optical 
disc, offset by about 30kms _1 from the main component. 
This might be the counterpart of the expanding ionised gas. 
Other than that, the neutral gas does not seem to follow the 
motions of the ionised gas. 

Altogether, the kinematic analysis confirms the sugges- 
tion we made in Sect. |3.2[ i.e., that we detected an expand- 
ing superbubble in the north-east of NGC 5408. However, 
the relatively low spectral resolution and the poor spatial 
coverage make it impossible to model this superbubble. 
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5 DISCUSSION 

Our analysis has shown that both galaxies show differences 
in the kinematic behaviour of the neutral and ionised gas 
components. Often, this can be interpreted as evidence for 
outflows of the ionised gas. Sometimes, these outflows are 
detected in both Ha and Hi. In this section, we discuss our 
results and make some predictions about the fate of the gas. 
We will also assess our results with respect to other studies 
of dwarf galaxies. Finally, we will come back to the peculiar 
H II region in the south of IC 4662. 



5.1 The outflows and the special case of IC 4662 

The outflow velocities vary from 20 up to 60kms -1 in both 
galaxies. In some cases we found an Hi counterpart for an 
expanding ionised gas structure, moving at about the same 
velocity. NGC 5408 shows little evidence for line-splitting in 
H I. The neutral gas distribution seems to be warped in the 
outer parts, which might indicate gas accretion, interaction 
with gas-rich companions or infall of H I clouds (see |Sancisi| 
et al.|2008 and references therein). However, NGC 5408 ap- 
pears to be quite isolated and no H I clouds could be found 
in its near vicinity on our deep Hi maps. Furthermore, gas 
infall would result in a lopsided gas distribution, i.e., more 
gas on one side than on the other, which is definitely not 
seen in NGC 5408. Therefore, the outflow scenario is much 
more likely. 

From the morphology and the kinematics, we suggest 
that NGC 5408 harbours at least one expanding superbub- 
ble. In comparison to IC4662, this galaxy shows outflows 
almost everywhere across the disc. 

For IC 4662 no XMM Newton observations are avail- 
able and the ROSAT X-ray All-Sky Survey does not show 
a source that could be associated with this galaxy. Never- 
theless, IC4662 is a very complex system: the comparison 
of the neutral and ionised gas kinematics revealed that it is 
very difficult to take the H I velocities as a reference value for 
the Ha velocities as the Hi velocity field is quite distorted. 
Especially the discrepancy in slit 3 (see Fig. 



12 



lower right 

panel) is probably due to the large-scale H I velocity field and 
not due to an outflow of the ionised gas. The same might be 
true for slit 2 (between 30" and 70") where the Hi veloci- 
ties again drop significantly. However, here, the ionised gas 
seems to follow that movement. If we just look at the veloc- 
ities of the ionised gas and try to define a velocity gradient 
across each slit, a velocity offset of about 20kms _1 blue- 
shifted can be seen, which coincides with the north-eastern 
filamentary area. From the morphology alone expanding gas 
is expected, but because of the distorted Hi velocity field, 
and the lack of spatial coverage and spectral resolution of 
the long-slit spectra, it is difficult to quantify the outflow 
velocities. 

As can be seen very nicely in NGC 5408 (Fig. 



13 lower 



row), the kinematics of the ionised gas are dominated by 
stellar feedback causing offsets from the rotation velocity of 
the galaxy. This implies that the Ha velocities can not be 
assumed to reflect the rotation of a galaxy; instead, we rely 
on the H I velocities. However, in the case of IC 4662, the H I 
velocity field is very distorted as well, which is why we are 
very careful about our detections. 

Looking at the large-scale structure, IC 4662 shows line- 



splitting of the H I line at many positions across the galaxy 
(Fig. [7]). As already mentioned above (Sect. 4.1), the two 
components are often separated by about 70kms _1 . Given 
the size of the optical galaxy in comparison to the distri- 
bution of the neutral gas, this is probably not caused by 
an interaction of stars and the ISM. As we see global line- 
splitting, the infall of single H I clouds is also highly unlikely. 
Instead, the line-splitting is probably a consequence of a 
large scale distortion like, e.g., a merger process. As IC4662 
is a very isolated galaxy, we may rule out any kind of cur- 
rent tidal interaction with another object. We think that it is 
more likely that two dwarf galaxies collided in the past and 
are now in the process of merging. This hypothesis explains 
the distorted kinematics we observe in IC4662, the two per- 
pendicular systems as observed in the Hi distribution, and 
also the Hi tail we found in the eastern part of the galaxy 
(see Fig. |3}. There is increasing evidence that interactions 
and mergers of dwarf galaxies trigger star-formation activ- 
ity in dwarf objects, but such features can only be detected 
when optical and radio observations are carefully compared 
L6pez-Sanchez et al.||2010 ). 



(e.g. 



5.2 Large FWHMs in IC 4662 and NGC 5408 

The Ha line profiles of both galaxies show very large 
FWHMs (see Fig.[l4|. In the case of IC4662, the FWHMs 
of the low intensity profiles are at about 70kms _1 . Even for 
high intensity profiles, this number does not decrease, which 
means that the high FWHMs are not an effect of low S/N. 
70 km s _1 are far too high to be explained by thermal broad- 
ening alone as this would require an electron temperature of 
125000°! Taking a reasonable value of the electron temper- 
ature of 10000 K gives a contribution of the thermal compo- 
nent of about 23kms _1 . This still leaves us with 66kms _1 
that have to be explained by other mechanisms. As shown 
in Sect. |4.4T7] the residuals after subtracting a single Gaus- 
sian from the observed line profile usually indicate that the 
Ha line is a superposition of at least two components. As 
the line has a symmetric shape, this may indicate expanding 
bubbles. 

NGC 5408 also shows quite high FWHMs. Here, it be- 
comes obvious that this result is also not an effect of 
the spectral resolution as the line widths remain equally 
broad when measured at higher spectral resolution (from the 
archival spectrum). In comparison to IC4662, the FWHMs 
are more scattered, although they are preferably located in 
a band centred at 70kms _1 , which is similar to the aver- 
age value found for IC4662. Additionally, the ionised gas 
observed in slit lb reaches FWHMs of up to lOOkms -1 . 
These high values have mainly been found in the western 
H II region, which is where the X-ray source is located. 

Detailed studies of the velocity dispersion in Hll re- 



gions have been carried out by Munoz-Tunon et al. (1996) 



and later by Martmez-Delgado et al. ( 2007 ) . Their diagrams 



show very similar patterns compared to ours: a narrow hor- 
izontal band and inclined bands at higher velocity disper- 
sion, which they interpret e as coming from faint expanding 
bubbles and shells. Indeed, we found the largest FWHMs 
to coincide with areas of filamentary structure (IC 4662) or 
potential superbubbles (NGC 5408). 
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5.3 Outflow or galactic wind? 

As we could not derive a reasonable rotation curve from the 
distorted velocity field of IC4662, we focus this analysis on 
NGC 5408. 

Following the method described in | van Eymeren et al.| 



(2009a), we compared the expansion velocities with the es- 



cape velocity of NGC 5408 in order to make predictions 
about the fate of the outflowing gas structures. The es- 
cape velocity was estimated by using the cored pseudo- 
isothermal halo model, which represents the shape of the 
rotation curves best, especially in the inner few kpcs (e.g., 
van Eymeren et al. 2009), which is where we detected the 



outflows. Only the circular velocity and the maximum halo 
radius r max have to be known in order to calculate the es- 
cape velocity. The circular velocity was estimated from the 
plateau of the H I rotation curve. We derived the escape ve- 
locity for two different radii. First, we took the outer radius 
of the Hi intensity distribution (9.5 kpc) as a lower limit for 
the halo size. As a second estimate we calculated the es- 
cape velocity using twice the Hi radius (19 kpc), which we 
consider to be a more physical assumption. The resulting 



curves are shown in Fig. 15 (dotted and solid lines respec- 



tively). The observed rotation curve is indicated by small 
grey triangles; the outflows with their distances from the 
dynamic centre and their velocities compared to the H I ve- 
locities are marked by large black triangles. A higher r max 
results in an increase of the escape velocity. As the radius 
of the dark matter halo r max is probably much higher than 
assumed, the escape velocities will be even higher than the 
ones we derived. The curve with r max = rm^ax gives there- 
fore a lower limit for the escape velocity values. As can be 
seen, the expansion velocities of the outflows, although par- 
tially quite significant, stay far below the escape velocity, 
even for the lower limit case. This is mainly due to the fact 
that they were found close to the dynamic centre where the 
escape velocities are high. 



5.4 Comparison with previous work 

We now compare our results for IC 4662 and NGC 5408 with 
similar studies from the literature. Many dwarf galaxies re- 
veal filaments and shell-like structures up to kpc-size on deep 
Ha images (see, e.g., |Hunter et aL 1993 |Bomans et al. 



1997: van Eymeren et al.|2007 ). Spectral line data show Ha 
line splitting, indicating outflowing gas, which is often as- 
sociated with these filaments and shells (e.g., Martin|19 98). 
Partly, the neutral gas is visibly affected by the expand- 
ing ionised gas as can, e.g., be seen in the starburst dwarf 
galaxy NGC 5253 flKobulnicky fe SkillmanpOO^ . The neu- 
tral gas is pushed by the ionised gas, which leads to the 
observed expansion velocities. On the other hand, the sur- 
rounding neutral gas decelerates the expanding bubble so 
that it might eventually be stalled by the neutral gas en- 
velope ( jKobulnicky fe S killman 2008). Figure [4] shows very 
nicely how the expanding gas can affect the neutral g as, as, 
e.g., NGC 5408 has a local H 1 minimum in the centre of the 
potential superbubble. A very similar morphology has been 
found in NGC 5253 ([Lopez- Sanchez eT"aLl[2010| . NGC 2366 



van Eymeren et al. 2009a) and NGC 4861 (van Eymeren 
et al.|2009b) also show H I minima coinciding with expand- 
ing ionised gas structures. 
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Figure 15. Outflow or galactic wind? The expansion velocities of 
the outflowing gas structures in NGC 5408 are compared with the 
escape velocities calculated by using the pseudo-isothermal halo 
model. The escape velocity is calculated for an isothermal halo 
of r max = 9.5 kpc (dotted line), which fits the size of the Hi dis- 
tribution, and r max = 19 kpc (solid line). The observed rotation 
curve is indicated by small grey triangles. The error bars repre- 
sent receding and approaching side. The expanding gas structures 
are marked by black large triangles. 



The expansion velocities of the galaxies in our sam- 
ple (including NGC 2366 and NGC 4861) vary between 
20kms _1 and 60kms _1 . These values are in good agree- 
ment with previous studies: [Heckman et al.| ( [T995| ) and |WestT 
moquette et al. ( 2008 ) found expansion velocities betwee n 50 



and lOOkrns" 1 in NGC 1569; Schwartz & Martin 



(2004), us- 



ing absorption line studies, found the expansion velocities of 
a sample of dwarf galaxies to be about 30kms _1 . An excep- 
tion seems to be the Wolf-Rayet blue compact dwarf galaxy 
He 2-10, where expansion velocities up to 250kms _1 have 
been detected ( jMendez et al.|1999 ). Note that more massive 
galaxies like M 82 sho w outflows with expansion velocities 
higher than 100 kms -1 (|Schwartz & Martin|2004||Strickland 



&; Heckman|2007 Westmoquette et al.|2009| ). In our sample, 
the highest velocities have been found in NGC 5408. 

How high do the expansion velocities have to be to turn 
an outflow into a galactic wind? Coming back to our sample 
(including NGC 2366 and NGC 4861, but excluding IC4662 
whose velocity field is too distorted to derive a rotation curve 
and the associated escape velocity), it could be shown that 
the expansion velocities of the outflows only account for 30 
to 50% of the escape velocities. Note that these velocities 
were calculated under the assumption that the radius of 
the dark matter halo r max equals the Hi radius. As r max 
is probably much higher, the resulting escape velocities will 
be higher, which therefore further decreases the possibility 
of a galactic wind. This result is consistent with studies of 
other dwarf galaxies (e.g., Martin|1998 Schwartz & Martin 
2004). An exceptio n is NGC 1569 which rotates very slowly 



with only 30 km s _1 ( Heckman et al. 1995 ) , but harbou rs out- 
flows with expansion velocities of up to 100 km s -1 ( |Heck-| 
man et al.|199 5; West moquette et al.|2008 ). Still, a galactic 
wind could not be confirmed, although the chances are much 
higher than in other dwarf galaxies. Another exception is 
He 2-10 (Men dez et al.||1999 ), which however, as mentioned 
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before, shows very fast expanding gas with velocities up to 
250 km s" 1 . 

At first glance, all of these non-detections of a galactic 
wind look like a contradiction with theoretical models like, 



e.g., the superbubble model by Mac Low & Ferrara (1999) 



Note, however, that outflowing gas is mostly detected close 
to the dynamic centre, where the escape velocities are very 
high. Only a few kpcs away from the dynamic centre, the es- 
cape velocity drops significantly (see, e.g., Fig. 15), increas- 



ing the chance of a galactic wind. As the surface brightness 
goes with n 2 e with n e being the electron density, a less dense 
gas will be much more difficult to detect. This implies that 
the outer regions of an outflow will easily escape the de- 
tection in emission line studies (e.g., Bomans et al. 2007), 
which means that we consistently miss the fastest moving 
gas. A study of the interstellar Nal absorption line (which is 
as an absorption line not affected by the n 2 e problem) does 
not help since Nal traces much cooler, denser gas, which is 
probably not as much accelerated as the diffuse HII. 

Additionally, we have to take into account the dynamic 
masses of our sample galaxies, which are of the order of a few 
times 10 9 M©, i.e., close to the upper limit of the mass range 
of dwarf galaxies. As Fig. 2 in Mac Low &; Ferrara (1999) 
shows, at total masses of 10 9 Mq, most of the gas is still 
bound to the disk when the energy input stops. This strong 
dependence of blow-out on the mass of the galaxy has al- 
ready been noted by |De Young &; Heckman| ( |1994[ ). Further- 
more, the central star clusters in, e.g., NGC 2366 ( |Drissen 
et al.|200ip or NGC 4861 (e.g., |Barth et al.|1994| ) are only a 
few Myrs old. If we assume that the current star formation 
episode is the driver of the outflows, it makes a superbubble 
blowout even more unlikely. 

It is also possible that the expansion velocity of the out- 
flows will never get close to or beyond the escape velocity 
due to the dark matter halo that slows down the gas that 
was expelled from the disc ( |Silich &; Tenor io-Tag le 1998). 
Their 2d calculations of multi-supernova remnants evolving 
in dwarf galaxies show that the total retention of the out- 
flowing gas is very likely for galaxies with an ISM mass of 
about 10 9 M . 

Altogether, this means that our observations (including 
NGC 2366 and NGC 4861) are in good agreement with other 
observations and with the simulations. All sample galaxies 
harbour filaments and shells that could be shown to expand, 
but no evidence for a galactic wind could be found. 

As the expansion velocities in IC4662 are comparable 
to those detected in the other three galaxies, we suppose 
that we would find similar results for this galaxy. However, 
if this galaxy is a merging system, processes with much more 
impact than star formation activity probably dominate the 
kinematic behaviour. 



5.5 The southern H n region in IC 4662 

The southern Hll region in IC4662 has been and still is 
a subject of great interest. It appears to be detached from 



the main body. Furthermore, Hidalgo- Gamez et al. (2001) 
pointed out that it differs significantly from the main body 
in oxygen abundance. Therefore, the question arises whether 
it is an own system interacting with IC4662 (e.g., H. Lee, 
private communication) or just an unusually placed OB as- 
sociation. Our deep Ha image shows that this bright Hll 



region is at least connected to the main complex by a chain 
of small Hn regions and diffuse filamentary gas structures 
(see Sect. [3] and Noeske et al. |2003| . Moreover, together 
with the main body this region is embedded into the high- 
est column density Hi. Additionally, its kinematics do not 
differ from the main body, but follow the rotation of the 
whole optical system. All these detections suggest that the 
southern Hll region is part of IC4662. 

Similar features have been found in the blue compact 
dwarf galaxy NGC 1140: a chain of Hll regions that em- 
anates south of the main body, and a core which is kinemat- 



ically decoupled from the rest of the galaxy (Hunter et al. 
1994; We stmoquette et al.|2 010). Therefore, this galaxy was 
suggested to be a merger remnant ( Hunter et al.||1994 ). A 
study of 20 starburst galaxies by |L6pez-Sanchez (2010) has 
shown that in fact almost all of them show signs of inter- 
action or of merging with small companions. They conclude 
that interactions play an important role in the triggering 
mechanism of the strong star- format ion activity observed in 
dwarf galaxies. 

Taking into account all the observational peculiarities 
and this similarity to NGC 1140, we definitely need to dis- 
cuss the detached H 11 region in IC 4662 in the context of a 
possible merger. If IC 4662 indeed consists of two merging or 
merged systems, the difference in the chemical composition 
might indicate that the main Hll region complex belongs 
to one galaxy and the southern Hll region to the second 
galaxy. Detailed mult i- wavelength observations of IC4662 
(and other blue compact dwarf galaxies) are on the way 
(see also Lopez-Sanchez et al.||2010 ). 



6 SUMMARY 

We analysed H 1 and optical data of the two peculiar dwarf 
galaxies IC 4662 and NGC 5408 to compare the morphology 
and the kinematics of the neutral and ionised gas compo- 
nents. 

In both galaxies, we detected outflows in Ha, sometimes 
with an Hi counterpart, which show expansion velocities of 
20 to 60kms _1 . The FWHMs of the Ha lines are as high as 
70 km s _1 , independent of the intensity of the lines. This sug- 
gests that other mechanisms apart from thermal broadening 
are taking place. In many cases, the Ha line is probably a 
superposition of at least two components, which indicates 
that expanding superbubbles might cause the high FWHMs. 
Higher spectral resolution data are necessary to study the 
outflows in more detail. 

We compared the expansion velocities measured in 
NGC 5408 with the escape velocity of the galaxy and con- 
clude that it is very unlikely that the gas will be blown away. 
This is in agreement with many previous studies. 

As IC4662 has a very distorted Hi velocity field, we 
did not perform an analysis of the fate of the gas. Never- 
theless, we detected outflowing ionised gas. The neutral gas 
often consists of two components, separated by 70kms -1 . 
We suggest that this galaxy is actually a merger of two inde- 
pendent dwarf systems. The most peculiar feature in IC 4662 
is the southern H 11 region, whose origin is still not clear. Al- 
though our analysis indicates that it is probably a part of 
the galaxy, it may also be the remnant of one of the two 
merged systems. 
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APPENDIX A: INTENSITY- WEIGHTED MEAN 
VS. HERMITE VELOCITY FIELD 

Here, we present the intensity-weighted mean (iwm) and the 



Hermite velocity field of IC 4662 (Fig. Al left and middle 
panels). By applying Gauss-Hermite h3 polynomials to a 
set of spectra, the peak velocities can be defined more accu- 
rately as this method is not biased by low S/N or blended 
lines. However, in the case of IC 4662 (and also in the case 
of NGC 5408), the H I profiles are either symmetric and sim- 
ple or the components can well be differentiated. This is 
confirmed if we subtract the Hermite velocity field from the 



iwm velocity field (Fig. Al right panel): the two velocity 
fields only differ in areas of very low S/N. The large scale 
velocity field remains the same. 
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Figure Al. A comparison of the intensity- weighted mean (iwm) velocity field (left panel) vs. the Hermite h3 velocity field (middle 
panel). The residuals after subtracting the Hermite from the iwm velocity field is shown on the right panel. 



